Quantum chaos algorithms and dissipative decoherence with quantum trajectories 
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Using the methods of quantum trajectories we investigate the effects of dissipative decoherence in 
a quantum computer algorithm simulating dynamics in various regimes of quantum chaos including 
dynamical localization, quantum ergodic regime and quasi- integr able motion. As an example we use 
the quantum sawtooth algorithm which can be implemented in a polynomial number of quantum 
gates. It is shown that the fidelity of quantum computation decays exponentially with time and 
that the decay rate is proportional to the number of qubits, number of quantum gates and per gate 
dissipation rate induced by external decoherence. In the limit of strong dissipation the quantum 
algorithm generates a quantum attract or which may have complex or simple structure. We also 
compare the effects of dissipative decoherence with the effects of static imperfections. 
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The main fundamental obstacles in realization of quan- 
tum computer are external decoherence and inter- 
nal imperfections. The decoherence is produced by cou- 
plings between the quantum computer and the external 
world (see e.g. review [2j). The internal imperfections 
appear due to static one-qubit energy shifts and resid- 
ual couplings between qubits which exist inside the iso- 
lated quantum computer. These imperfections may lead 
to emergence of quantum chaos and melting of quan- 
tum computer eigenstates 0, • The effects of unitary 
quantum errors produced by decoherence and imperfec- 
tions on the accuracy of quantum algorithms have been 
studied by different groups using numerical modeling 
of quantum computers performing quantum algorithms 
with about 10 — 20 qubits. The noisy errors in quantum 
gates produced by external decoherence are analyzed in 

IB, IE H H E3, EH E3, EE E3, EE 13 while the errors 

induced by internal static imperfections are considered in 
E 13 El El EJ El E3 . The analytical treatment UK 
based on the random matrix theory allows to compare the 
accuracy bounds for these two types of errors for quan- 
tum algorithms simulating complex quantum dynamics. 

In fact, a convenient frame for investigation of quan- 
tum errors effects in quantum computations is given by 
models of quantum chaos 24]. Such models describe a 
quantum dynamics which is chaotic in the classical limit 
and which has a number of nontrivial properties including 
dynamical localization of chaos, quantum ergodicity and 
mixing in phase space (see e.g. |24|)- ^ nas been shown 
that for many of such models the quantum computers 
with n q qubits can simulate the quantum evolution of an 
exponentially large state (e.g. with N = 2 Uq levels) in a 
polynomial number of elementary quantum gates n g (e.g. 
with rig = 0{n 2 q ) or n g = 0(nj*)). The quantum algo- 
rithms are now available for the quantum baker map l2al . 
the kicked rotator the quantum sawtooth UtI Il9| 
and tent 21] maps, the kicked wavelet rotator |18| . the 
quantum double- well map 10]. Their further generaliza- 
tion and development gave quantum algorithms for the 
Anderson metal-insulator transition |20(, electrons on a 



lattice in a magnetic field and the kicked Harper model 
[2^ . The quantum algorithm for the quantum baker map 
has been implemented experimentally with a NMR based 
quantum computer [27j . 

However till now the quantum chaos algorithms have 
been used only for investigation of unitary errors effects. 
This is always true for internal static imperfections but 
the external decoherence generally leads also to dissipa- 
tive errors. The first step in the analysis of dissipative de- 
coherence in quantum algorithms has been done in |2£| on 
a relatively simple example of entanglement teleportation 
along a quantum register (chain of qubits). After that 
this approach has been applied to study the fidelit y d ecay 
in the quantum baker map algorithm [2£| . In [23, ^ ne 
decoherence is investigated in the Markovian assumption 
using the master equation for the density matrix written 
in the Lindblad form lU- Already with Tin = 10 — 20 
qubits in the Hilbert space of size N = 2™ the numeri- 
cal solution of the exact master equation becomes enor- 
mously complicated due to a large number of variables in 
the density matrix which is equal to N 2 . Therefore the 
only possibility for numerical studies at large n q is to use 
the method of quantum trajectories for which the num- 
ber of variables is reduced to N with additional averaging 
over many trajectories. This quantum Monte Carlo type 
method appeared as a result of investigations of open 
dissipative quantum systems mainly within the field of 
quantum optics but also in the quantum measurement 
theory (see the original works 0, EE EE 133). More re- 
cent developments in this field can be find in Hi, SHI). 

In this paper we investigate the effects of dissipative 
decoherence on the accuracy of the quantum sawtooth 
map algorithm. The system Hamiltonian of the exact 
map reads [X7J Uil 



H 8 (n, 0) = Th 2 /2 + kV(0) s (t ~ m ) 



(1) 



Here the first term describes free particle rotation on 
a ring while the second term gives kicks periodic in 
time and n = —id/d0. The kick potential is V{0) = 



2 



-(0 - tt) 2 /2 for < < 2tt. It is periodically re- 
peated for all other so that the wave function ip(9) sat- 
isfies the periodic boundary condition ip(0) = ij){6 + 2tt). 
The classical limit corresponds to T — > 0, /c — > oo with 
if = fcT = const. In these notations the Planck con- 
stant is assumed to be ft = 1 while T plays the role of an 
effective dimensionless Planck constant. 

The classical dynamics is described by a symplectic 
area-preserving map 

n = n + k(6 -tt), = Q + Tn. (2) 

Using the rescaled momentum variable p = Tn it is easy 
to see that the dynamics depends only on the chaos pa- 
rameter K = kT. The motion is stable for —4 < K < 
and completely chaotic for K < —4 and K > (see 
and Refs. therein). The map (0) can be studied on the 
cylinder (p G (— oo, +oo)), which can also be closed to 
form a torus of length 27rL, where L is an integer. 

The quantum propagation on one map iteration is de- 
scribed by a unitary operator U acting on the wave func- 
tion ip: 

^=U^j = U T U k ^ = e -^/2 e ^kV(0)^ (3) 

The quantum evolution is considered on N quantum mo- 
mentum levels. For N = 2 nq this evolution can be im- 
plemented on a quantum computer with n q qubits. The 
quantum algorithm described in performs one itera- 
tion of the quantum map JSJ in n g — 3n^ + n g elementary 
quantum gates. It essentially uses the quantum Fourier 
transform which allows to go from momentum to phase 
representation in n q (n q + l)/2 gates. The rotation of 
quantum phases in each representation is performed in 
approximately n 2 q gates. Here we consider the case of 
one classical cell (torus with L = 1 when T = 2tt/N) [l^ 
and the case of dynamical localization with N levels on 
a torus and K ~ 1, k = K/T ~ 1 19]. Here and below 
the time t is measured in number of map iterations. 

To study the effects of dissipative decoherence on the 
accuracy of the quantum sawtooth algorithm we follow 
the approach with the amplitude damping channelused 
in [2^. The evolution of the density operator p(t) of 
open system under weak Markovian noises is given by 
the master equation with Lindblad operators L m (m = 
l,--- ,n q ): 

P=-\ i H effP ~ P H lf f\+J2 L mP^m , (4) 

m 

where the system Hamiltonian H s is related to the ef- 
fective Hamiltonian H e ff = H s — ih/2 s ^ Jm L\ n L rn and 
m marks the qubit number. In this paper we assume 
that the system is coupled to the environment through 
an amplitude damping channel with L m = am^/Y, where 
a m is the destruction operator for m— th qubit and the 
dimensionless rate T gives the decay rate for each qubit 
per one quantum gate. The rate T is the same for all 
qubits. 




FIG. 1: Probability distribution W n over momentum eigen- 
states n in the quantum sawtooth map Q at time t = 30. The 
quantum evolution is simulated by the quantum algorithm 
with n q = 6 qubits in presence of dissipative decoherence. 
The dissipation rate per gate is F — 0.001 and the map pa- 
rameters are k = \/3, K — y/2 with the total number of states 
N = 2 nq = 64. The full curve represents the exact solution of 
the Lindblad equation 0). Symbols show the result of quan- 
tum trajectories computation with the number of trajectories 
M = 20(+), M = 50(o), M = 200(x), M = 1000(A). The 
initial state is n — 0. The logarithm is natural. 

This evolution of p can be efficiently simulated by av- 
eraging over the M quantum trajectories which evolve 
according to the following stochastic differential equation 
for states \ifj a ) (a = 1, • • • , M)\ 

\dr) = -iH s \r)dt + \ J2((LlLmh (5) 

m 

-LiL m )\r)dt + ( / Lm - 1 ) \^ a ) dN m , 

m (LtnLm)^ J 

where represents an expectation value on l^") and 
dN m are stochastic differential variables defined in the 
same way as in 29] (see Eq.(10) there). The above 
equation can be solved numerically by the quantum 
Monte Carlo (MC) methods by letting the state \^ a ) 
jump to one of Lm,\^ a ) /\Lm, \^ a ) I states with prob- 
ability dp m = \L m \ilj a )\ 2 dt [2<j or evolve to (1 — 
iH eff dt/K)\^ a )/^/l - ^2 m dp m with probability 1 - 
^2m dpm- Then, the density matrix can be approximately 
expressed as 

1 M 

P (t) « (\m)W)\) M = mI1 infirm , (6) 

where ()m represents an ensemble average over M quan- 
tum trajectories \ip a (t)). Hence, an expectation value of 
an operator O is given by (O) = Tr(Op) w (O)m- 




FIG. 2: (Color online) First top row: classical phase space 
distribution obtained from the classical sawtooth map with 
a Gaussian averaging over a quantum cell (N — 256 quan- 
tum cells inside whole classical area; see text). Second row: 
the corresponding Husimi function for the quantum sawtooth 
map at n q = 8 and r = 0. Third row: the Husimi function ob- 
tained with M = 50 quantum trajectories in presence of dissi- 
pative decoherence with rate T = 0.0005 and n q — 8. Fourth 
bottom row: same as for the third row but with n q = 10. Here 
K = -0.5, T = 2tt/N corresponding to L = 1 and N = 2 n< ? 
quantum states in the whole classical area. Columns show 
distributions averaged in the time intervals: < t < 9 (left), 
40 < t < 49 (middle), 90 < t < 99 (right). The initial state is 
n w O.liV. Color represents the density from blue/black (0) 
to red/gray (maximal value). 

For the quantum sawtooth algorithm the dissipative 
noise is introduced in the quantum trajectory context 
(Eq. (0) after each elementary quantum gate and cal- 
culated by the MC methods. The same physical process 
can be described by density matrix theory. The evolu- 
tion of density matrix after single iteration of quantum 
sawtooth map is described by 

p' = U k U T pU f T Ul . (7) 

To include the dissipative noise effects the density matrix 
further evolves according to Eq. (|3J with H s = between 
consecutive quantum gates composing Uk and Ut- 
To test the accuracy of the method of quantum trajec- 
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FIG. 3: Fidelity / as a function of iteration time t. The upper 
two curves are for T = 0.0005 and the lower two curves are 
for T = 0.001. Here M = 50, n g = 8, k = 2 n *K/2it and 
K = —0.5 (full curves) K = 0.5 (dotted curves), respectively. 
The initial state is \n = 0). The inset shows the fidelity decay 
rate 7 as a function of T e ff = n q n g Y with n q — 4, 6, 8. Here 
K = 0.5 (+) and K = -0.5(A), respectively. The straight 
line is the best fit 7 = 0.08 r e //. 



tories we compare its results with the exact solution of 
the Lindblad equation for the density matrix p (gj). The 
comparison is done for the case of dynamical localiza- 
tion of quantum chaos and is shown in Fig. 1. It shows 
that the dynamical localization is preserved at relatively 
weak dissipation rate Y. It also shows that the quan- 
tum trajectories method correctly reproduces the exact 
solution of the Lindblad equation and that it is sufficient 
to use M = 50 trajectories to reproduce correctly the 
phenomenon of dynamical localization in presence of dis- 
sipative decoherence of qubits. 

To analyze the effects of dissipation rate T in a more 
quantitative way we start from the quasi-integrable case 
K = —0.5 with one classical cell L = 1 (T = 2tt/7V). The 
classical phase space distribution, averaged over a time 
interval and a Gaussian distribution over a quantum cell 
with an effective Planck constant, is shown in Fig. 2 in the 
first top row (there are TV = 2 Uq quantum cells inside the 
whole classical phase space). Such Gaussian averaging 
of the classical distribution gives the result which is very 
close to the Husimi function in the corresponding quan- 
tum case at T = (Fig. 2, second row). We remind that 
the Husimi function is obtained by a Gaussian averaging 
of the Wigner function over a quantum H cell (see |38j for 
details). In our case the Husimi function h(6,n) is com- 
puted through the wave function of each quantum trajec- 
tory and after that it is averaged over all M trajectories. 
The effect of dissipative decoherence with T = 0.0005 is 
shown in the third row of Fig. 2. At T = the phase 
space distribution remains approximately stationary in 
time while for F > it starts to spread so that at large 
times the typical structure of the classical phase space 
becomes completely washed out. This destructive pro- 
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FIG. 4: (Color online) Dependence of the IPR £ on time t 
at T = 0.001 and M = 50 shown by full curves for n q = 4 
(green/gray curve), n g = 6 (blue/black curve), n q — 8 (black 
curve), bottom to top respectively. The dashed curves show 
the same cases at Y = (bottom dashed curve is for n q = 4, 
top dashed curve is for n q = 6 and n g = 8 where £ values 
are practically identical). Here the initial state is \n = 0) and 
k = V3, K = y/2. 

cess becomes more rapid with the increase of the number 
of qubits even if T remains fixed (Fig. 2, fourth bottom 
row). One of the reasons is that T is denned as a rate 
per gate and the number of gates n g = 3n 2 + n q grows 
with n q . However, this is not the only reason as it shows 
the analysis of the fidelity decay. 

The fidelity / of quantum algorithm in presence of dis- 
sipative decoherence is defined as 

f{t) = (Mt)\p(t)\Mt)) * ^ E \(Mt)W(t))\ 2 , (8) 

a 

where |^o(^)) is the wave function given by the exact 
algorithm and p(t) is the density matrix of the quantum 
computer in presence of decoherence, both are taken after 
t map iterations. Here, p(t) is expressed approximately 
through the sum over quantum trajectories (see also 29]). 

The dependence of fidelity f(t) on time t is shown in 
Fig. At relatively short time t < 50 the decay is ap- 
proximately exponential f(t) ~ exp(— 7^). The decay 
rate 7 is described by the relation 

j = CT eff = Cn q n g T , (9) 

where C = 0.08 is a numerical constant (see Fig. |2l inset). 
The important result of Fig. 01 is that the decay of f(t) 
is not very sensitive to the map parameters. Indeed, it is 
not affected by a change of K which significantly modify 
the classical dynamics which is quasi-integrable at K = 
—0.5 and fully chaotic at K = 0.5. Another important 
result is that up to a numerical constant the relation Q 
follows the dependence found in 29] for the quantum 




FIG. 5: (Color online) Dependence of the ratio of £ to its 
value £0 in the ideal algorithm on the dissipative decoherence 
rate F for n q = 4 (green/gray o), n q = 6 (blue/black x), 
n q = 8(black +) from bottom to top. The values £ and £0 are 
averaged in the time interval 30 < t < 40. Other parameters 
are as in Fig. 0] 

baker map. This shows that the dependence © is really 
universal. Its physical origin is rather simple. After one 
gate the probability of a qubit to stay in upper state 
drops by a factor exp(— T) for each qubit (we remind that 
T is defined as a per gate decay rate). The wave function 
of the total system is given by a product of wave functions 
of individual qubits that leads to the fidelity drop by a 
factor exp(— Cn q T) after one gate and exp(— Cn q n g Y) 
after n g gates leading to the relation (0. In principle, 
one may expect that the decay of f(t) is sensitive to a 
number of qubit up states in a given wave function since 
there is no decay for qubit down states. However, in 
a context of a concrete algorithm this number varies in 
time and only its average value contributes to the global 
fidelity decay. 

The result © gives the time scale tf of reliable quan- 
tum computation in presence of dissipative decoherence. 
On this scale the fidelity should be close to unity (e.g. 
/ = 0.9) that gives 

t f ^l/(n q n g T); N g = 1/ (n q T) . (10) 

Here N g = n g t is the total number of quantum gates 
which can be performed with high fidelity (/ > 0.9) at 
given n q and T. The comparison with the results ob- 
tained for static imperfections [2]| of strength e shows 
that for them N g drops more rapidly with n q : N g ~ 
l/(e 2 n q n g ). Therefore the static imperfections destroy 
the accuracy of quantum computation in a more rapid 
way compared to dissipative decoherence. 

It is also interesting to analyze the effects of dissipative 
decoherence on the dynamical localization. For that, in 
addition to the probability distribution as in Fig. ^ it 
is convenient to use the inverse participation ratio (IPR) 
defined as £ = l/EJV'nl 4 - V E„ KKMVl 2 where 
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FIG. 6: (Color online) Each panel shows the Husimi distribu- 
tion for the quantum sawtooth map algorithm with K = 1, 
T — 2n/2 nq and n q = 8. The top three rows show the cases 
with the rate F = 0.01, T = 0.05, and F = 0.1, respectively 
from top to third row. The initial state is \n = 60) and 
M — 50. The distribution is averaged in the time interval 
< t < 9 (left column), 40 < t < 49 (middle column), 
90 < t < 99 (right column). The fourth bottom row shows 
the distribution for another initial state \n = 0) averaged 
in the time interval 90 < t < 99 for F = 0.01 (left panel), 
r = 0.05 (middle panel) , and r = 0.1 (right panel) (compare 
with the right column of top three rows). Color represents 
the density from blue/black (0) to red/gray (maximal value). 

(| • • • |)m notes the average over M quantum trajectories. 
This quantity is often used in the problems with localized 
wave functions. In fact £ gives an effective number of 
states over which the total probability is distributed. The 
dependence of £ on time t is shown in Fig. 2] It shows 
that in presence of dissipative decoherence the dynamical 
localization is destroyed. Indeed, at large t the value 
of £ grows with n q while for the ideal algorithm it is 
independent of n q . The physical meaning of this effect 
is rather clear. As in Fig. [3| the dissipative decoherence 



introduces some noise which destroys localization. 

However, there is also another effect which becomes 
visible at relatively large F. It is shown in Fig. which 
gives the ratio of £ to its value £o in the ideal algorithm. 
Thus, at small F the ration £/£o grows with the increase 
of T while it stars to drop at large F. This is a manifesta- 
tion of the fact that in absence of algorithm the dissipa- 
tion drives the quantum register to the state \n = 0) with 
all qubits in down state. Even in presence of the quan- 
tum algorithm this dissipative effect becomes dominant 
at large F leading to a decrease of the ratio £/£o- 

The dissipative effect of decoherence at large values of 
T is also clearly seen in the case of quantum chaos er- 
godic in one classical cell (L = 1). At large F the Husimi 
distribution relaxes to the stationary state induced by 
dissipation \n = 0) (third row in Fig. at F = 0.1). In a 
sense this corresponds to a simple attractor in the phase 
space. The stationary state becomes more complicated 
with a decrease of F (second row in Fig. [Slat F = 0.05). 
And at even smaller F = 0.01 the stationary state shows 
a complex structure in the phase space (top first row in 
Fig.|SJ). It is important to stress that this structure is in- 
dependent of the initial state (bottom row in Fig.|SJ). In 
this sense we may say that in such a case the dissipative 
decoherence leads to appearance of a quantum strange 
attractor in the quantum algorithm. Of course, this sta- 
tionary quantum attractor state is very different from the 
Husimi distribution generated by the ideal quantum al- 
gorithm. However, it may be of certain interest to use 
the dissipative decoherence in quantum algorithms for 
investigation of quantum strange attractors which have 
been discussed in the context of quantum chaos and dis- 
sipation (see e.g. Refs. [33, E3, El| ) . At the same time 
we should note that the dissipation induced by decoher- 
ence acts during each gate that makes its effect rather 
nontrivial due to change of representations in the map 

In conclusion, our studies determine the fidelity de- 
cay law in presence of dissipative decoherence which is 
in agreement with the results obtained in |29| for a very 
different quantum algorithm. This confirms the univer- 
sal nature of the established fidelity decay law. These 
studies also show that at moderate strength the dissipa- 
tive decoherence destroys dynamical localization while a 
strong dissipation leads to localization and appearance 
complex or simple attractor. The effects of dissipative 
decoherence are compared with the effects of static im- 
perfections and it is shown that in absence of quantum 
error correction the later give more restrictions on the 
accuracy of quantum computations with a large number 
of qubits. 

This work was supported in part by the EU IST-FET 
project EDIQIP. 
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